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Thermodynamic and thermoeconomic investigations of a small-scale integrated gasification solid oxide 
fuel cell (SOFC) and Stirling engine for combined heat and power (CHP) with a net electric capacity of 
120 kW e have been performed. Woodchips are used as gasification feedstock to produce syngas, which is 
then utilized to feed the anode side of the SOFC stacks. A thermal efficiency of 0.424 LHV (lower heating 
value) for the plant is found to use 89.4 kg/h of feedstock to produce the above mentioned electricity. 
Thermoeconomic analysis shows that the production price of electricity is 0.1204 $/kWh. Furthermore, 
hot water is considered as a by-product, and the cost of hot water is found to be 0.0214 S/kWh. When 
compared to other renewable systems of similar scales, this result shows that if both SOFC and Stirling 
engine technology enter the commercialization phase, then they can deliver electricity at a cost that is 
competitive with the corresponding renewable systems of the same size. 

© 2014 Elsevier Ltd. All rights reserved. 


1. Introduction 

Due to the ever-increasing demand for more efficient power 
production and distribution, the main areas of electricity produc¬ 
tion research and development are efficiency enhancements and 
pollutant reduction, especially carbon dioxide reduction. There is 
an increased interest in developing a distributed system of smaller- 
scale facilities rather than using one large-scale capacity at a 
specified location, which would allow electricity and heat to be 
produced and distributed close to the end user, minimizing the 
costs associated with transportation [1,2], 

SOFC stacks will soon enter the commercialization phase, while 
small Stirling engines are approaching this phase. Therefore, it 
would be interesting to combine these two technologies into a 
single system, quantifying the benefits of each system to establish a 
new technology. Together with an integrated gasification plant that 
gasifies woodchips in a two-step gasification process, electricity and 
heat could then be produced in an environmentally friendly way. 

SOFCs are one of the most promising types of fuel cells, partic¬ 
ularly for energy production. They are expected to produce clean 
electrical energy at high convention rates, with low noise and low 
pollution emissions [3], 

Due to the high exhaust temperatures of SOFCs caused by the 
high operating temperature of the cells and the fact that the fuel 
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utilization in the fuel cell never reaches 100% [4], the unreacted fuel 
can be combusted in a burner, which in turn produces off-gases 
with even higher temperature that are perfectly suited for use in 
a heat engine, such as a Stirling engine, for the production of power 
and heat for domestic purposes. 

The numerous studies on SOFC-based power systems in the 
open literature have suggested high thermal efficiencies. However, 
the majority of these studies use gas turbines as the bottoming 
cycle (see, e.g., Refs. [5—7]). A steam turbine has also been used as a 
bottoming cycle [8], resulting in high plant efficiency. Only a few 
studies have been carried out with a Stirling engine as the bot¬ 
toming cycle when a fuel cell cycle is used as topping cycle (see, e.g., 
Refs.! 1,2,9,10]). At present, using Brayton and Rankine cycles as the 
bottoming cycles seems to be more applicable because of the 
maturity of these technologies. As development trends suggest that 
the SOFC operating temperature will decrease in the future, using a 
gas turbine as bottoming cycle will be less beneficial. 

Integrated gasification SOFC systems have also been studied for 
some time (see, e.g., Refs. [11—15]). However, integrated biomass 
gasification SOFC—Stirling CHP (combined heat and power) plants 
have not been investigated in the open literature, forming the basis 
of this study. 

The present work is an analytical study that includes both 
thermodynamic and thermoeconomic investigations of systems 
with integrated gasification of woodchips, where the syngas is used 
as fuel for a SOFC plant, which also functions as the topping cycle 
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for a Stirling engine by contributing heat from the exhausted off¬ 
gases, see Fig. 1. The system net capacity is 120 l<W e , which is 
suitable for decentralized CPH plants. The gasifier model used for 
the analysis is based on a downscaled version of the Viking two- 
stage gasifier build at DTU-MEK and now located at DTU Riso. 
This gasifier produces a clean syngas that can be directly fed into an 
SOFC (for more information on the gasifier plant, see Refs. [16—18]). 
The SOFC is based on theoretical equations with empirical co¬ 
efficients from an experimental arrangement, while the Stirling 
engine parameters are chosen by fitting a validated feasible engine 
in terms of construction. 

2. Methodology 

The thermodynamic results in this study are obtained from the 
Dynamic Network Analysis (DNA) simulation tool (see, e.g., 
Ref. [19]). This software is a result of an ongoing development 
process in the Thermal Energy Section of the Mechanical Depart¬ 
ment of the Technical University of Denmark, which began with a 
master’s thesis project [20], Since this time, the program has been 
continuously developed to be generally applicable, covering unique 
features and hence supplementing other simulation programs. 

The program includes a component library, thermodynamic 
state models for fluids and standard numerical solvers for differ¬ 
ential and algebraic equation systems. The component library 
models include heat exchangers, burners, turbo machinery, dryers, 
decanters, energy storages, valves and controllers, among others. 
The thermodynamic state models for fluids cover most of the basic 
fluids and such compounds as ash and tar for use in energy system 
analyses. 

DNA is a component-based simulation tool, meaning that the 
model is formulated by connecting components with nodes and 
adding operating conditions to build up a system. Next, the physical 
model is converted into a set of mathematical equations and solved 
numerically. The equations include mass and energy conservation 
for all components, and the nodes represent the relationships 
among the thermodynamic properties of the fluids in the system. 
The total mass balance and energy balance for the entire system is 
also included to account for heat loss and heat exchange between 
different components. In addition, the components include a 
number of constitutive equations representing their physical 
properties, e.g., heat transfer coefficients for heat exchangers and 
isentropic efficiencies for compressors and turbines. The program is 
written in FORTRAN, and users may also contribute additional 
components and thermodynamic state models to the libraries. 
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Fig. 1. Block scheme of the plant. 


For the thermoeconomy, the general theory is combined with 
the governing component cost equations to create a system of 
linear equations, which is then solved in Engineering Equation 
Solver (EES). 


2.1 Modeling of the SOFC stacks 


The SOFC model developed in this investigation is based on the 
planar type developed by DTU Riso and Tops 0 e Fuel Cell 
(topsoefuelcell.com). The model was validated against experi¬ 
mental data in the range of 650—800 °C (the range of operating 
temperatures), as described in Ref. [21]. For clarity, this model is 
described briefly as follows. The model is assumed to be zero¬ 
dimensional, thus enabling the calculation of complicated energy 
systems. In such modeling, one must distinguish between elec¬ 
trochemical modeling, the calculation of cell irreversibility (cell 
voltage efficiency) and the calculation of the species compositions 
at the outlet. For electrochemical modeling, the operational voltage 
(£ C eii) is found to be 

^cell = ^Nernst — AE a ct — AE 0 hm — A£ C onc, (1) 

where ENernst. AE a ct. A£ 0 hm and A£ C onc are the Nernst ideal revers¬ 
ible voltage, activation polarization, ohmic polarization and con¬ 
centration polarization, respectively. Assuming that only hydrogen 
is electrochemically converted, then the Nernst equation can be 
written as 


^Nernst 


~ Ag f I RT ]n ( P»2.tOtVP^ \ 
n e F + n e F \ p H2 o )' 


( 2 ) 


PH 2 ,tot = Ph 2 + Pco + 4Pch 4 ! (3) 

where Ag? is the Gibbs free energy (for H 2 reaction) at standard 
pressure. Because the electrochemical oxidation rate of H 2 is much 
higher than that of CO, then it can be assumed that hydrogen is the 
only species that can convert electrochemically (see Refs. [22,23]). 
In the above equations, p Hl and p H2 o are the partial pressures for H 2 
and H 2 O, respectively. 

The activation polarization can be evaluated from the Butler— 
Volmer equation, which is isolated from other polarizations to 
determine the charge transfer coefficients and exchange current 
density from the experiment using the curve fitting technique (see, 
e.g., Refs. [24,25]). 


AEact — AE ac t c + A£, 


2RT 

~T 


sinh 


act,a 

1 (i + in 


2 i. 


0 ,a 


+ sinh 


1 / I + In 


2 i, 


0 ,c 


(4) 


The ohmic polarization depends on the electrical conductivity of 
the electrodes as well as the ionic conductivity of the electrolyte 
(see, e.g., Refs. [26,27]). This polarization is also calibrated against 
experimental data for a cell with anode, electrolyte and cathode 
thicknesses of 600 |im, 50 pm and 10 pm, respectively. 

The concentration polarization is dominant at high current 
densities for anode-supported SOFCs, wherein insufficient amounts 
of reactants are transported to the electrodes and the voltage is 
then reduced significantly. Again, the concentration polarization is 
calibrated against experimental data by introducing the anode 
limiting current (see, e.g., Refs. [28,29]), in which the anode 
porosity and tortuosity are also included, among other parameters. 

The fuel composition at the anode outlet is calculated using the 
Gibbs minimization method, as described in Ref. [30], Equilibrium 
at the anode outlet temperature and pressure is assumed for the 
following species: H 2 , CO, CO 2 , H 2 O, CH 4 and N 2 . Thus, the Gibbs 
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minimization method calculates the compositions of these species 
at the outlet by minimizing their Cibbs energies. The equilibrium 
assumption is fair because the methane content in this study is very 
low. The methane content is mainly depended on the kinetic pa¬ 
rameters rather than the chemical equilibrium and its reaction rate 
is fast [31], 

To calculate the voltage efficiency of the SOFC cells, the power 
production from the SOFC (Psofc) depends on the amount of 
chemical energy fed to the anode, the reversible efficiency (ijrev). 
the voltage efficiency (r; v ) and the fuel utilization factor (L/ F ). It is 
defined in mathematical form as 

Psofc = (LHV H2 ii H2 in + LHV C0 h C0in + LHV CH4 h CH4 in )?] rev 7j v Lf F , 

(5) 

where U F is a set value and r\ v is defined as 


working fluid flows back through the regenerator, and the regen¬ 
erator rejects heat to the working fluid. The temperature of the 
working fluid increases in the second isochoric process (4 —> 1). 

The model of the Stirling engine provides a rather conservative 
engine in terms of efficiency, among other properties. It is mainly 
used to show the applicability of the technology as a functioning 
bottoming cycle for the usage of a hot exit flow, which in this case is 
generated by combustion after achieving a rather high temperature 
in the SOFC topping cycle. The model parameter inputs are selected 
such that construction is feasible, e.g., infinite surface areas of heat 
exchangers are not allowed. The heat source used in the analysis is 
the combustion product gases from the catalytic gas burner, while 
water used for domestic purposes is used as the sink. 

In this study, a pseudo-Stirling cycle, which provides better 
agreement with engine performance data, is adopted [33], The 
power output for the Stirling engine is modeled as 


Vv = 


AT+ell 

pNernst 


( 6 ) 


The reversible efficiency is the maximum possible efficiency, 
defined as the relationship between the maximum electrical en¬ 
ergy available (change in Gibbs free energy) and the fuel’s lower 
heating value (LHV) as follows (see, e.g., Ref. [32]): 
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LHV fuel 


(7) 


(Ag f ) 


fuel 


gf 


h 2 o 

gf 


- ISf) -jia 


co 2 




gf 


co 2 


2 gf 


h 2 o 


- gf 


yH 2 ,in 
yco,i 
~2 (I f 


CH, 


Tea,, in- 

( 8 ) 

The partial pressures are assumed to be the average of the inlet 
and outlet temperatures, 


Pj = ( Ss ^ !Us )p j= {H 2 ,C0,CH 4 ,C0 2 ,H 2 0,N 2 }, 


Po 2 


/y 0 2 ,out +y 0 2 ,in 

l- 2 - 


Pc- 


(9) 


In addition, equations for the conservation of mass (with molar 
flows), conservation of energy and conservation of momentum are 
also included in the model. The reliability of the models presented 
here is justified in Ref. [2] with different fuels such natural gas, 
ethanol, and methanol. 


^Stirling = Ppcy (Qhigh ~ Qlo 

where Qi oss is defined as 
Qloss = Qhigh(l “ Pmec.stirl) i 


( 10 ) 


( 11 ) 


where 7) me c,stiri is the mechanical efficiency of the Stirling engine 
and Qhigh is the amount of heat the Stirling engine absorbs from the 
heat source. The polytropic efficiency 7) pcy is defined as [33] 


Vpcy 


(1 - RV 1 ~ T ) - £(RV T_1 - 1) 

(1 - RV 1_T ) + (1 - ?)(1 - estiri). 


( 12 ) 


where RV and f s tiri are the reversibility factor for the Stirling engine 
and the effectiveness of the internal heat exchangers in the engine, 
respectively. The constant y is equal to 1.667, and 'C is defined as the 
temperature of the cooler gas over the heater gas, 


? = 


* cooler,gas 


+ 273.15 


^heater,; 


gas ' 


273.15 


where Theater gas is 
Theater,gas = Theater, wall — ATjugp, 

T'cooler.wall = T'water, inlet + 0.66667(AT wa ter), 

T'cooler,gas = 7 C ooler.wall + AT| 0W . 


(13) 

(14) 

(15) 

(16) 


2.3. Modeling of the gasifier 


2.2. Modeling of the Stirling engine 

The Stirling engine is noted for its quiet operation and the ease 
with which it can make use of almost any heat source. Stirling 
engines are referred to as external combustion heat engines and 
operate based on a regenerative closed power cycle using helium, 
nitrogen, air or hydrogen as the working fluid. An ideal regenerative 
Stirling cycle consists of four processes in one cycle. First, the 
working fluid absorbs the heat from a high-temperature reservoir 
and undergoes isothermal expansion (process 1 —> 2). Second, the 
hot working fluid flows through a regenerator, and the regenerator 
absorbs heat from the hot working fluid. Thus, the temperature of 
the working fluid decreases in an isochoric process (2 —» 3). Third, 
the working fluid rejects heat to a low-temperature reservoir and 
undergoes isothermal compression (3 —> 4). Lastly, the cold 


The Viking two-stage gasifier is modeled and used for the 
analysis in this investigation. The Viking gasifier plant is a 75 l<W t h 
biomass gasifier using an autothermal (air-blown) fixed bed 
gasifier. This gasifier is built and developed by the Biomass Gasifi¬ 
cation Group at the Technical University of Denmark [16]. Wood 
pellets are used as feedstock. The pellets are first dried to remove 
surface moisture and pyrolyzed in the first reactor, after which the 
pyrolytic gases (600 °C) are fed into an autothermal downdraft 
gasification reactor. The produced exhaust gases after the gasifier 
are used to heat the reactor for the drying and pyrolysis processes 
(see Fig. 3). Partial oxidation creates a high-temperature tar¬ 
cracking zone in between the pyrolysis and the char gasification 
step (zone) by injecting preheated air. Char is gasified in the fixed 
bed, and H 2 0 and C0 2 are the gasifying agents in the char gasifi¬ 
cation reactions. The gasifier operates at approximately atmo¬ 
spheric pressure. 
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The gasifier is modeled by a simple Gibbs reactor, where the Gibbs 
free energy is at a minimum when chemical equilibrium is achieved. 
This characteristic is used to calculate the gas composition at a spe¬ 
cific temperature and pressure without taking the reaction paths into 
account [30] and will be briefly explained below. The Gibbs free 
energy of a gas (assumed to be a mixture of k perfect gases) is 

C=X> \s° + RT ln(n,p)j, (17) 

i=t 

where g°, R and T are the specific Gibbs free energy, universal gas 
constant and gas temperature, respectively. Furthermore, each 
atomic element in the inlet gas is in balance with the outlet gas 
composition, which shows that the flow of each atom must be 
conserved. For N elements, this balance is expressed as 

k w 

= zE "m,out A mj forj=l,JV, (18) 

i = 1 m = 1 

where the N elements correspond to H2, O2, N2, CO, NO, CO2, steam, 
NH3, H2S, SO2, CH4, C, NO2, HCN (hydrogen cyanide), COS (carbonyl 
sulfide), Ar and ash (Si02) in the gasifying process. A mj is the 
number of atoms of element j (H, C, O or N) in each molecule of the 
entering compound i (H 2 , CH 4 , CO, CO2, H2O, 0 2 , N 2 or Ar), while A,j 
is the number of atoms of element j in each molecule of the exiting 
compound m (H2, O2, N2, CO, NO, CO2, steam, NH3, H2S, SO2, CH4, C, 
NO2, HCN (hydrogen cyanide), COS, Ar or ash). The minimization of 
the Gibbs free energy is formulated by introducing a Lagrange 
multiplier, fi, for each of the N constraints obtained. The expression 
can be simplified as 


Heat exchangers are all assumed to counter the flow, and it is 
assumed that all of the energy is transferred from one side to the 
other, neglecting heat losses. Both the logarithmic mean tempera¬ 
ture difference (LMTD) and the effectiveness-number of transferred 
units (e-NTU) methods are used to calculate the flows in the heat 
exchangers, depending on the type [34], 

To model the methanator, which is used to increase the methane 
content in the fuel, a methanation process is used. This process is 
mainly expressed by the exothermic reaction of CO and H 2 , which 
are reformed into CH 4 and steam under high pressure. However, 
some other minor reactions will also occur. 


CO + 3H 2 —> CH 4 + H 2 0. (23) 

In the catalytic gas burner, unused fuel is reformed in a highly 
exothermic process, which is modeled based on the general equa¬ 
tion for the burning of hydrocarbons in oxygen, 

C n H m + (n + ™)o 2 ^nC0 2 + (™)h 2 0, (24) 


where n and m denote the amount of the hydrocarbon present in 
the reaction. Furthermore, for all reforming components, the 
analysis software uses the Gibbs free energy for the different 
compounds. Thus, chemical equilibrium is obtained by minimizing 
the Gibbs free energy for the reactions. 

The pump's power consumption is obtained by volumetric and 
pressure states as 


Wpump 


muin (Pout e; Pin)' 


V 


pump 


(25) 


J = i 


— Gtot,out H~ ^ ^ Mj ( 'y ^ ^i,outAj y ^ ^m,inA 


mj 


i = 1 


m = 1 


where m, p, v and 77 are the mass flow, pressure, specific volume 

(19) 

(m 3 /kg) and efficiency of the pump, respectively. 


Furthermore, by setting the partial derivation of this equation 
with respect to n,- out to zero, the function can be minimized as 


clt/l _ dtjtot.cn 


+ E HjAy = 0 for i = 1 , k. 


i =i 


( 20 ) 


>g- out + RT ln ( n i,outPout) + E PjAj = 0 for ' = 1 ’ k - 
i= i 


Thus, a set of k equations is defined for each chemical compound 
leaving the system. 

In the previous study, Rokni [15] compared the computed dry 
wood gas composition with the measured dry wood gas composi¬ 
tion from Viking gasifier and showed that the computed results 
agreed very well with the measured data. 


2.4. Modeling of the other components 


The compressor power consumption is modeled based on the 
definition of the isentropic efficiency rjisentropic and mechanical ef¬ 
ficiency TJmechanical as 


Pisntropic 


'W _ c±h' 0 

W ~ M o 


V — T 
1 exit 1 m 


Tp \it — Fu¬ 


ll! J compressor 


( 21 ) 


2.5. Cost model and cost of components 

To create the cost model, general theory from Ref. [35] is used to 
obtain the costs of each separate stream in the energy system and 
to determine the cost of the produced electricity and domestic hot 
water, which is considered a by-product. Here by-product means 
that it is a secondary product which has no influence on financial 
results and therefore they do not receive allocation of joint cots. The 
system is modeled using different components, denoted k, with 
different associated costs. The cost rates of each component are 

h ■ ^CI -OM 

Cpk — C Fk -\-Z k Z k . (2b) 

where C Pk is the total cost rate associated with the product of the 
component, while C Fk is the cost rate associated with the fuel 

(inflow) to the component. Z k ' and Z° M are the cost rates associated 
with the annual contribution to investment and operation and 
maintenance, respectively, for component k divided by the annual 
operation time of the system in hours. Thus, to determine the cost 
of the entire system, all component cost balances need to be taken 
into consideration: 

Cp = C F+ zZ+Z™. (27) 


^mechanical — 


Ah 0 


W 


compressor 


( 22 ) 


where the mark ' represents the compression process without 
isentropic losses and W, h and T represent work, enthalpy and 
temperature. 


The cost rates of the product and fuel are derived from the 
exergy flow E and the specific exergy cost c and are provided below 
for the exergy transfer rate £), £ e . power, and heat transfer, 
expressed in $/kWh. 

Q = CjE; = c^rhjei), C e = c e E e = Ci(m e ee), (28) 
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Cw — Cq — CqEq. 


(29) 


Table 1 

Estimates of rates for total investment cost. 


The component cost rates are determined using the total capital 
investment (TCI) cost, which is found by defining the purchase cost 
of each component (PEC). The PEC cost equations are obtained from 
Ref. [ 36 ] and are in $. 

/sOFC = /stack,SOFC + /aux.SOFC > 

/stack,SOFC = ( W ce // 7rD cell i cell)( 2 -96T ce | 1 - 1907), (30) 

/aux.SOFC = 0.10/soFC, 


where cylindrical cells with a diameter of 0.005 m and a length of 
0.12 m are assumed. For the SOFC inverter, the cost is 


/sOFC, inverter — 10 


We, 

500 


0.7 


(31) 


Total capital investment 

Percentage of PEC 

Direct costs 


(a) Onsite costs 


Purchased equipment installation 

33% 

Piping 

35% 

Instrumentation and controls 

12% 

Electrical equipment and materials 

13% 

(b) Offsite costs 


Civil, structural and architectural work 

21% 

Service facilities 

35% 

Indirect costs 


Engineering and supervision 

8% 

Construction and contractors’ profit 

15% 


(W e | is in units of W, and the other W s, e.g., compressors, are in 
units of kW). 

The cost of the Stirling engine is assumed from Ref. [1 ] with the 
assumption that the total TCI is 2200 $/kW installed. Furthermore, 
this price represents the present cost and is therefore expected to 
decrease when Stirling engines reach commercialization levels. 
The compressor purchase cost is obtained from 
/ ■ \ 0 67 

/compressor = 91,562 f ^ j . (32) 

The purchase cost of the heat exchangers is obtained from 

^ '“(ns®)" 8 ' (33) 

and the heat exchanger area is calculate by 


The purchase costs of the splitters, mixers and valves are 
neglected due to the scaling. Thus, the costs related to the other 
components are considered to be part of the direct costs related to 
the investment, e.g., piping, instrumentation and controls. 

The TCI of PEC for each component k is 



149\ 

+ TooJ 


23 

+ Too 



(38) 


where the relationships among the direct and indirect costs related 
to the investment used are shown in Table 1. 

From the TCI of each component, the two cost rates Z k and Z k 
can be derived. Z k is expressed as 


7 CI - 
2- k — 


h 


h/year 


(39) 


2*hex 


mill 


(34) 


where the TCI of component k is amortized in n years by the an¬ 
nuity factor/as 


where the values for k are 35 W/m 2 1< for gas-gas heat exchangers 
and 135 l<W/m 2 K for gas—liquid heat exchangers. A7i m is the log¬ 
arithmic mean temperature difference of the heat exchangers. 

For the purchase cost of the fuel reformer, (the gas cleaner and 
methanator), the following equation is used 

/pR = 13 °(ow) + 3240(V , p R ) 0,4 + 2128 V P r, (35) 


/ x TCI, 


/ 


( q<" +cp >-l ^ 

W-Dsrv 




-i -i 


where qi is the interest factor obtained from 


<?i = 




(40) 


(41) 


where the finned area is assumed to be 10% of the volume of the 
fuel reformer. V/r is obtained from the inlet mass flow and the 
reaction time of the process. The reaction time is assumed from 
Ref. [ 37 ] to be 1 s for the gas cleaner and 0.5 s for the methanator. 
The gas burner is assumed to be a catalytic burner, and its cost is 
therefore assumed to be the same as that for the fuel pre-reformer. 

The pump purchase cost is obtained from 


/pump — 442/w pll mp) 1.41/,, 


/; — 1 + 


1-0.8 A 

1 pump / 


(36) 


The gasifier purchase cost is assumed from Ref. [ 38 ] as follows, 
considering the small scale of the gasifier: 

/ . n 0.67 

/gasifier = 1600 (ihclrywood J i ( 37 ) 


where the dry wood mass flow is given in kg/h. 


CP, int and r, are the construction period, interest rate and rate of 
inflation, respectively, and are provided in Table 2. For simplifica¬ 
tion, the lifetime of the whole system is assumed to be 20 years. The 
construction period is assumed to cover the time from the initial 
planning until the final completion and commissioning. The in¬ 
terest rate and inflation rates are assumed from Ref. [ 39 ], with 
Denmark used as the area of reference. 


Table 2 

Economical parameters assumed for the cost rates. 


Parameter 

Value 

Annual operating hours 

7500 [h/year] 

Interest rate 

3% 

Inflation rate 

2% 

Equipment lifetime 

20 [years] 

Construction period 

1 [yearj 

Maintenance factor 

1.1 
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Furthermore, with M as the maintenance factor, Z k can be 
obtained from 


■yOM 


= Z k x M. 


(42) 


This gives the combined contribution to the cost rates associated 
with the investment of the components. 

Since thermoeconomic equations for most components 
(compressor, turbine, heat exchanger, etc.) are well-known 35], 
only equations for the gasifier, SOFC and Stirling engine are 
described here in detail. 

A fuel cell integrated with a bottoming cycle can be described by 
different cost models, varying by product and fuel allocation [35,36], 
In this study the exergy difference between the outgoing used fuel 
and the inlet reformed gas is considered as fuel. Electric power and 
flue gas are considered as product. Cost and exergy balance together 
with auxiliary equations for SOFC (referring to Fig. 2a) can be found 
in Table 3. 

Similarly, for the Stirling engine (referring to Fig. 2b), exergy and 
cost balances together with auxiliary equations can also be found in 

Table 3. 

For the gasifier (referring to Fig. 2c), only the syngas is allocated 
as the product. The fuel is made up of steam-air mixture and dried 
woodchips. Cost and exergy balances are expressed in Table 3. In 
order to make the system determined, the auxiliary (steam 99) sets 
the cost of ash disposal in €/l<Wh equal to zero. 

The average price for woodchips is assumed to be approximately 
45 DKK/GJ [40], which is then converted into $/kWh using the 
currency rate conversion of 1 $ — 5.7DKK. Transportation costs for 
the wood are neglected and it is assumed to be obtained locally and 
therefore have low associated costs. Transportation cost is almost 
insignificant compared to other cost. In many countries wood are 
transported through rivers. 

By combining all of the component equations, a linear system is 
created in EES to obtain the costs of all streams of the system. 
Furthermore, for all components, a single fuel and product are 
defined for each component so that the cost balances for all com¬ 
ponents could be found. No heat loss is considered in the thermo¬ 
dynamic analysis; hence, E q (defined in Eq. (29)) is neglected. Thus, 
the system could be solved by adding a number of auxiliary equations 


Table 3 

Energy and cost balance for the gasifier, SOFC and Stirling engine. 


Gasifier 


SOFC 


£74 +E 2 = £3 + E 99 + Ed, gasifier 
C74E74 + C2E2 +Z gasifler = C3E3 + C99E99 

El, gasifier = Egg 
C99 = 0 

£21 - £22 = E28 - E27 + Eel + Eo.SOFC 

C2I E2I - C22E22 +Z S 0 FC = c 2 P,E' 2 H - C27E27 + Cp e | S0FC P e | 

El, sofc = 0 


C 21 = c 22 

/- — ^28^28 ~ C27E27 

Pel, SOFC E 28 -E 2 7 

Stirling engine £ 40 _ £ 4I = f 43 - £ 42 + P el + £ D stirling 

C40E40 - C41E41 + ^Stirling = c 43^48 - ^42^42 + c Pel,StirlingPel 
^L,Stirling = 0 
c 40 = C41 

C _ . _ C42E42—C43E43 

L Pei, Stirling 


based on the assumptions that the fuel flowing through a component 
can be used in a later stage, which means that the fuel had the same 
unit cost at the inlet and outlet. If the product of a component is 
composed by two or more streams, then the unit cost of those 
streams are assumed to be the same. Thus, for the heat exchangers 
and dryer, the specific exergy cost of the hot streams is considered to 
be equal. Furthermore, for compressors and pumps, the specific 
exergy cost of the inlet ambient air or ambient water is assumed to be 
zero. For the gas cleaner and gasifier, the disposal costs of sulfur and 
ash are neglected, meaning that the specific exergy cost is zero. The 
splitters’ outlet streams are assumed to be equal. For SOFC, the as¬ 
sumptions are 1) that the exergy difference of the inlet fuel and the 
outlet fuel is the fuel consumption of the SOFC and hence had equal 
specific exergy costs and 2) that the flue gas is considered as waste, 
and hence the specific exergy cost is zero. In this case, the fuel for 
Stirling bottoming cycles still had costs related to the streams, and if 
more fuel is left in the stream (due to the lower fuel utilization factor 
of the SOFC), then this assumption is reasonable, and in this way, the 
bottoming cycle is not considered a full regenerative cycle. For the 
Stirling engine, two auxiliary equations are needed. First, the outlet 
stream on the hot side is considered as waste and hence zero. Addi¬ 
tionally, the water heater utilizing this waste is assumed to recover 


Fuel 


Air 


21 


27 


Anode 


SOFC 


Cathode 


Off-fuel 


22 

►El. 


28 


Off-air 



a) 


Dried woodchips 


Air and Steam 


74 


Gasifier 


99 | Ash 

c) 


b) 


Wood gas 


Fig. 2. (a) SOFC, (b) Stirling engine and (c) gasifier streams schemes. 
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energy from the waste. The second assumption is that the products of 
the Stirling components are mechanical power from the engine, and 
the heat absorbed by the heat source thus leads to 


c P,mec 


c cold,in^cold,in c cold,outsold,out 


^cold,in 


-E, 


cold,out 


(43) 


Electricity is considered the main product of the plant, and heat 
is considered the secondary product, which will be present as long 
as electricity is produced. It is thus assumed that the internal power 
consumption is covered by the production; hence, the internal 
power cost will be equal to the cost of the produced electricity: 


Price,,i = 


£p,SOFC c P,SOFC + £p,Stirling c P, Stirling 


£p,SOFC + £| 


P,Stirling 


(44) 


3. Plant configuration 

The system investigated herein is a small-scale CHP consisting of 
an integrated biomass gasification plant with an SOFC system as the 
topping cycle and a Stirling engine with a water heater as the 
bottoming cycle (see Fig. 3). 

Woodchips are fed into the system to produce syngas. This is 
performed in a two-step process, wherein the first step is the 
drying and pyrolysis of the feedstock and the second is fixed bed 
gasification, where the pyrolyzed feedstock is gasified by steam and 
air as gasification agents. Although Ref. [18 reported that the 
produced syngas is clean enough to be fed directly into the SOFC 
without additional fuel processing, a gas cleaner is introduced for 
the removal of H 2 S present in the syngas and bring its concentra¬ 
tion below 1 ppm level. The gas cleaner is assumed to work at a 
temperature of 250 °C. 

For the topping SOFC cycle, ambient air at 15 °C is compressed to 
the working pressure of the SOFC (normal pressure) before being 
heated in the cathode air preheater (CP) to the cathode inlet tem¬ 
perature of 600 °C. The cathode preheater uses some of the SOFC 
off-air for the heating. The off-air is split in two streams; one 
entering the CP and the other entering the catalytic burner (CB). For 
the anode side, the cleaned syngas is first pumped to compensate 


for the pressure drop along the way. Next, the syngas is sent to the 
methanator, wherein the molar fraction of CH 4 in the gas is 
increased from approximately 0.01 to nearly 0.05 at the expense of 
the molar fractions of H 2 , and CO and steam, while the molar 
fractions of N 2 and CO 2 also increase. This will not affect the SOFC’s 
electrical production in any particular way; however, because the 
reformation is highly exothermic, less heat needs to be extracted 
from the SOFC off-fuel to heat the incoming fuel to the SOFC fuel 
inlet temperature of 650 °C in the anode preheater (AP). The 
reformation allows the Stirling engine to receive a larger amount of 
heat to be used due to the higher temperature of the fuel when 
entering the CB, causing the combustion processes to occur at a 
higher temperature. The CB is implemented because not all of the 
fuel is reacted in the SOFC stacks due to fuel utilization. The 
entering temperatures are essential requirements for the proper 
functioning of the SOFC stacks in terms of both initiating the 
chemical reactions and avoiding cell thermal fractures. 

Second, a larger portion of CH 4 in SOFC causes endothermic 
internal reforming, which leads to the use of less air for cooling to 
maintain the SOFC operating temperature at 780 °C. Therefore, the 
workload of the cathode compressor/air blower will decrease. 

For the bottoming cycle, a Stirling engine is implemented. The 
Stirling engine utilizes the combustion products leaving the CB as 
the heat source. For the heat sink, water is used with an incoming 
temperature of 20 °C and an outgoing temperature of 60 °C so that 
it can be used, for example, as domestic hot water for room heating, 
at a temperature that is sufficient to address bacterial problems 
(e.g., Legionella [41]) and to be useful for domestic purposes. The 
heat remaining after the Stirling engine is used for domestic water 
heating. Water is constrained in the same manner as the heat sink, 
and the combustion products leaves the system into the environ¬ 
ment at approximately 95 °C, which is high enough to avoid 
corrosion problems. 


4. Results and discussion 

The main operating parameters for the plant are presented in 
Table 4. The ambient conditions are assumed to be 1 bar and 15 °C. 
The presented values are the final values used after the 
optimization. 



Fig. 3. System layout for the integrated gasification SOFC-Stirling-based CHP plant. 
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Table 4 

System operating input parameters. 


Parameter 

Value 

Woodchip temperature 

15 °C 

Pyrolysis and dry wood temperature 

150 °C 

Gasifier temperature 

800 °C 

Gasifier pressure 

1 bar 

Gasifier pressure drop 

0.005 bar 

Gasifier carbon conversion factor 

1 

Gasifier non-equilibrium methane 

0.01 

Steam blower isentropic efficiency 

0.8 

Steam blower mechanical efficiency 

0.98 

Steam temperature in steam loop 

150 °C 

Wood gas blower isentropic efficiency 

0.7 

Wood gas blower mechanical efficiency 

0.95 

Gas cleaner pressure drop 

0.0049 bar 

Cathode compressor air intake temperature 

25 °C 

Compressor isentropic efficiency 

0.7 

Compressor mechanical efficiency 

0.95 

SOFC operating temperature 

780 °C 

Anode inlet temperature 

650 °C 

Cathode inlet temperature 

600 °C 

Pressure drop, anode side 

0.02 bar 

Pressure drop, cathode side 

0.055 bar 

SOFC fuel utilization rate 

0.675 

Number of cells per stack 

74 

Number of stacks 

160 

Heat exchanger pressure drop 

0.01 bar 

Pinch temperature CP 

20 °C 

Burner ratio inlet outlet pressure 

0.97 

Stirling engine heater wall temperature 

600 °C 

Stirling engine AT high 

125 °C 

Stirling engine ATi ow 

60 °C 

Regenerator RV loss factor, Stirling engine 

1.44 

Heat exchanger efficiency, Stirling engine 

0.98 

Stirling engine mechanical efficiency 

0.8 

Stirling engine heat sink inlet temperature 

20 °C 

Stirling engine heat sink outlet temperature 

60 °C 

Water pump efficiency 

0.95 

Water heater inlet water temperature 

20 °C 

Water heater outlet water temperature 

60 °C 

Off gas temperature (after water heater) 

95 °C 


rather early stage. The SOFC has a power output of 98.8 l<W e , while 
the Stirling engine provides 26.9 l<W e of power. The internal power 
consumption is 5.8 kW, which is mainly due to the cathode air 
compressor. The high power consumption for the cathode air 
compressor is due to the cooling effect of this compressor, which is 
needed to maintain the SOFC temperature at the desired level. The 
reaction inside the SOFC is highly endothermic and therefore needs 
a relatively large flow of air to cool the cells. 

The thermal efficiency of the topping SOFC cycle is 0.329 LHV, 
which is somewhat low for a SOFC system (see Table 6). However, 
the entire plant has a thermal efficiency of 0.424 LHV, which is a 
decent value for such a small-scale system. The implementation of 
the bottoming Stirling engine gives a remarkable increase in plant 
efficiency of 28.9%. 

The Stirling engine heat sink produces nearly two thirds of the 
domestic water heating. The two water heaters together produce 
127.33 kW th . 

It is also assumed that the plant’s internal electrical consump¬ 
tion is covered by its production, meaning that the 120 l<W e pro¬ 
duction of electricity is the net production, 

fnet = PsOFC + ^Stirling ~ 2 T ei.n,consumedi (45) 

n 


where Psofc and Pstiriing are the power production from the SOFC 
and Stirling engine, respectively, and P e i,n,consumed is the consump¬ 
tion of power from the nth component. The thermal efficiency of 
the plant is calculated by the net power production of the SOFC and 
the Stirling engine compared to the fuel input: 


^th,plant 


P net 

difuel x LfWfuel 


(46) 


4.1. Thermodynamic investigation 


To achieve a 120 kW e output of electrical energy, the gasifier needs 
an input of approximately 89.4 kg/h, which leads to a syngas pro¬ 
duction of approximately 176.4 kg/h. Thus, this amount of biomass 
must be obtained from agriculture or a cultivation area and provided 
to the unit. The syngas molar fraction is presented in Table 5. 

Because ambient air is used as one of the gasifying agents, a 
large portion of unusable N 2 is present in the syngas. This issue is 
addressed by the application of a methanator prior to SOFC without 
using anode recirculation or external steam supplement. On the 
other hand, this leads to a larger mass flow, which is also beneficial 
for the Stirling engine operation. For the SOFC, however, a large 
amount of N 2 and steam causes a concentration polarization at a 


Table 5 

Syngas and methanated syngas molar fraction composition. 


Compound/molar fraction 

Wood 

Syngas 

Methanated syngas 

Hydrogen 

0.041 

0.253 

0.220 

Nitrogen 

0.001 

0.288 

0.310 

Carbon monoxide 


0.172 

0.086 

Carbon dioxide 


0.116 

0.186 

Water (liquid) 

0.332 



Water (gas) 


0.158 

0.147 

H 2 S (xlO 6 ) 


45.71 

- 

Methane 


0.010 

0.049 

Argon 


0.003 

0.004 

Oxygen 

0.293 



Carbon 

0.326 



Sulfur (xlO 6 ) 

133.6 



Ashes 

0.006 




The most important parameters for this investigation are the 
woodchip mass flow, number of SOFC stacks and SOFC cell utili¬ 
zation factor. The woodchip mass flow indicates the cultivation area 
to be allocated to provide the mass flow needed by the fuel. The 
number of SOFC stacks is directly related to the SOFC purchased 
cost and therefore investment cost, while the utilization factor af¬ 
fects the amount of off-fuel (fuel remaining after the SOFC stacks) 
that would be available for the bottoming cycle (the Stirling engine 
in this case). The lower utilization factor means that more fuel will 
be available for the Stirling engine; thus, the engine will produce 
more power. That is, the utilization factor affects the cooperation 
between the two cycles, namely, the SOFC plant and Stirling engine. 


Table 6 

Plant output for the initial operating parameters. 


Parameter 

System output 

Feedstock consumption 

89.4 kg/h 

Syngas produced 

176.4 kg/h 

Power output, SOFC 

98.8 kW 

Power output, Stirling engine 

26.9 kW 

Total power consumption 

5.8 kW 

Thermal efficiency, SOFC cycle 

0.329 

Thermal efficiency, plant 

0.424 

Percentage increase when adding Stirling cycle 

28.9% 

Component 

Produced heat 

Stirling engine 

53.5 

Water heater 

73.83 

Total hot water production 

127.33 
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a) 


b) 


Fig. 4. Plant thermal efficiency and net power production as function of woodchip mass flow with fuel utilization of 0.65 and 0.8: a) 100 stacks and b) 150 stacks. 


Following the discussion above, the first parameter to be 
investigated is the woodchip mass flow. To study the system per¬ 
formance with different woodchip mass flows, the initial analysis is 
performed on a system using 100 and 150 SOFC stacks with 74 cells 
per stack. The fuel utilization of the SOFC cells is initially assumed 
to be 0.8, see Ref. [4], The final utilization factor will be decided 
later. The plant efficiency and net power production as functions of 
woodchip mass flow are shown in Fig. 4 for two different utilization 
factors as 0.8 and 0.65. 

As shown, the thermal efficiency of the system tends to decrease 
as the mass flow in the system increases, and at a certain level of 
fuel mass flow, it tends to decrease rapidly. For the case of 100 
stacks, the power increases as the woodchip mass flow increases 
until reaching a maximum, after which it decreases. This behavior 
is also observed for the case with 150 stacks, but the woodchip 
mass flow is much higher. For other values of the SOFC utilization 
factor than 0.8, the results will be slightly different but the overall 
conclusions remain the same. 

Another important parameter is the SOFC utilization factor. To 
study this parameter, the woodchip mass flow is fixed to produce a 
net power of 120 kW at a utilization factor of approximately 0.7. 
Fig. 5 shows the variation of the plant efficiency and net power 
production when the SOFC utilization factor is changed. The results 
are shown for two different numbers of stacks (100 and 150) to 
study the plant performance. 

Several interesting points can be concluded from this figure. A 
utilization factor exists for which the plant efficiency is highest: 
0.625 for 100 stacks and 0.65 for 150 stacks. This maximum point is 


unique for a given number of stacks and must be calculated using 
the number of stacks chosen based on an economic analysis. 
Another issue to be mentioned is that the net power production 
and plant efficiency decrease sharply after a certain utilization 
factor. This point is found out to be approximately 0.8 for 100 stacks 
and 0.83 for 150 stacks. At this point, the concentration losses in the 
SOFC cells dominate, and as a result, the cell voltage decreases 
significantly. This point will obviously be shifted to the right when 
the number of stacks is increased. Increasing the number of stacks 
decreases the fuel amount for each stack when the total fuel mass 
flow remains constant. 

As mentioned above, the number of stacks is also an important 
issue to study. For this reason, several calculations are performed, 
and the results are illustrated in Fig. 6. For comparison purposes, 
two utilization factors are selected: 0.8, which is a rather high 
value, and 0.65, which is the optimum value when 150 stacks are 
used. Again, the fuel mass flow is fixed for these cases. For the case 
with a utilization factor of 0.8, increasing the stack number from 
100 to approximately 120 sharply increases the plant efficiency and 
power, and further increases in the stack numbers increase plant 
performance significantly. Eventually, the amount of fuel per stack 
becomes too large, and the ionic concentration reaches its limit. 
Increasing the stack number distributes the fixed fuel to more 
stacks and reduces the concentration to below its limit. 

For the case with a utilization factor of 0.65, increasing the stack 
numbers slightly increases plant efficiency and power. The problem 
of ionic concentration does not exist for such this lower utilization 
factor. Another conclusion from these results is that the plant 




a) 


b) 


Fig. 5. Plant thermal efficiency and net power production as function of utilization factor when the woodchip mass flow is fixed: a) 100 stacks and b) 150 stacks. 
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Fig. 6. Plant thermal efficiency and net power production as a function of the number of stacks when the woodchip mass flow is fixed: a) utilization factor = 0.8 and b) utilization 
factor = 0.65. 


efficiency for the case of a utilization factor of 0.65 is higher than 
that for the case of 0.8. However, the power produced with a uti¬ 
lization factor of 0.8 could be greater than that produced with a 
utilization factor of 0.65 if the number of stacks is high enough 
(more than approximately 160). The number of stacks is directly 
associated with the plant investment cost, and therefore a limit 
must be chosen to avoid high investment costs and thereby high 
prices for the produced electricity. The higher the number SOFC 
stacks, the higher the associated investment cost would be. By 
closer inspection of Fig. 6, choosing 160 SOFC stacks results in 
higher plant efficiency than choosing 120 stacks. On the other hand, 
the investment cost of 160 SOFC stacks is significantly lower than 
that of 200 stacks. Preliminary test results showed that the selec¬ 
tion of 160 stacks would be reasonable from the perspectives of 
plant efficiency and investment cost. For this reason, 160 stacks are 
used for the cost investigation, even though other values could be 
selected. 


4.2. Thermoeconomical investigation 

Based on the results from the thermodynamic analysis and the 
cost equations of each component along with the economic pa¬ 
rameters, the operation cost per hour [$/h] is obtained (see Table 7). 
The operation cost is denoted Z tot , which is the component in¬ 
vestment cost over its lifetime. The operation and maintenance 
costs are denoted Z 01 and Z° M , respectively. 


Table 7 

System and component cost rates [$/h]. 


Component 

Z 01 [$/h] 

Z° M [$/h] 

2“ [$/h] 

SOFC 

0.8344 

0.9179 

1.752 

Inverter 

0.0071 

0.0078 

0.015 

AP 

0.0301 

0.0332 

0.063 

Methanator 

0.0356 

0.0392 

0.074 

Syngas compressor 

0.0554 

0.0609 

0.116 

Gas cleaner 

0.0572 

0.0629 

0.120 

CP 

0.5794 

0.6374 

1.217 

Cathode compressor 

0.1032 

0.1136 

0.217 

Catalytic gas burner 

0.1030 

0.1134 

0.216 

Stirling engine 

0.3965 

0.4362 

0.833 

Water pump 

0.0033 

0.0036 

0.007 

Water heater 

0.0431 

0.0474 

0.090 

Gasifier 

0.7145 

0.7859 

1.500 

GAP 

0.0629 

0.0692 

0.132 

Steam generator 

0.0573 

0.0631 

0.120 

Steam loop blower 

0.0308 

0.0338 

0.065 

Total operating cost 

2.9495 

3.2444 

6.194 


As shown in the table, the largest cost is associated with the 
SOFC (1.752 $/h), here with 160 stacks and 74 cells in each stack, 
followed by the gasifier (1.5 $/h), the cathode air preheater (CP) 
(1.217 $/h), which is a relatively large heat exchanger, and then the 
Stirling engine (0.833 $/h). The costs are relatively low for the other 
components. The entire plant’s operating cost is found to be 
9.447 $/h with the current inputs and assumptions. 

The system’s sensitivity to the fuel price is illustrated in Fig. 7, for 
which the system costs presented in Table 7 are used. As shown, 
increasing the fuel price by 0.02 $/kWh increases the produced 
electricity cost by approximately 0.05 $/kWh. For the system with 
the current price of woodchips described in the cost model section, 
a production cost of 0.1215 $/kWh is found. The cost of the pro¬ 
duced domestic hot water (the cost of DHW) is found to be 
0.0219 $/kWh. This low cost is due to the assumption that the cost 
of DHW is the cost rate of the water heater and the exergy differ¬ 
ence of the “inflow” and “outflow” of the WH (water heater) and 
the Stirling engine. Hence, DHW is considered a byproduct of the 
plant because electricity is assumed to be the main product, and 
therefore the cost of the DHW will be the utilization of that heat. 

There are several options for reducing the cost of the system. As 
suggested in Table 7, the SOFC is by far the component that con¬ 
tributes most to the total operating cost. Because the number of 
stacks is a major contributor to the cost (see Eq. (38)), reducing the 
number of stacks will reduce the component investment cost 
significantly. However, as shown in Fig. 6, reducing the number of 
stacks will reduce the thermal efficiency of the plant, and because 
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Number of stacks 

Fig. 8. Cost of produced electricity over the number of stacks in the SOFC. 


the fuel cost is relatively low, the produced electricity cost will 
decrease with as the number of stacks decreases, see Fig. 8. 

As illustrated in Fig. 8, the produced electricity cost is minimized 
when the number of SOFC stacks is approximately 120. This finding 
suggests that with the current pricing, it would be more profitable 
to reduce the number of stacks to approximately 120 and increase 
the amount of feedstock to the gasifier, which in turn decreases the 
plant thermal efficiency but also the cost of electricity. Thus, the 
number of SOFC stacks can be determined by thermoeconomical 
optimization. Decreasing the number of stacks to 120 results in a 
different plant performance, which is shown in Table 8. 

The gasifier must be fed with 91.33 kg/h of wood to achieve a 
net power production of 120 kW. This amount of feedstock will 
lead to a syngas production of 180.05 kg/h, which means that the 
air intake into the gasifier is nearly equal to that of wood for the 
gasification processes. After the methanation of the syngas (to 
slightly increase the methane content), the power production of 
the SOFC plant will be approximately 98.48 kW e . As explained 
earlier, the methanator increases the amount of methane in the 
fuel, which in turn results in endothermic internal reforming. As a 
result, a smaller amount of air in the cathode side will be needed 
to cool the SOFC to maintain the SOFC operation temperature. 
Thus, the power consumption of the compressor will also be 
reduced. The total internal power consumption will be approxi¬ 
mately 6.10 kW. The power production from the Stirling engine is 
approximately 27.62 l<W e , leading to a plant thermal efficiency of 
0.415 (LHV-based). The integrated gasification and SOFC topping 
cycle will have a thermal efficiency of 0.32 (LHV-based). Thus, 
applying the Stirling engine as the bottoming cycle increases the 
plant thermal efficiency by nearly 28%, which is significant. This 


Table S 

Plant performance by thermoeconomical optimization. 


Parameter 

Value 

Fuel consumption 

91.33 kg/h 

Syngas produced 

180.05 kg/h 

Power output, SOFC 

98.48 kW 

Power output, Stirling engine 

27.62 kW 

Total internal power consumption 

6.10 kW 

Thermal efficiency, SOFC cycle 

0.3195 

Thermal efficiency, plant 

0.4149 

Increase when adding bottoming cycle 

29.88% 

Produced heat 

130.52 kW 

Cost of electricity 

0.1204 S/kWh 

Cost of DHW 

0.0214 S/kWh 

Total TCI over installed capacity 

3432.97 $/kW 


increase is due to the low thermal efficiency of the integrated 
gasification and SOFC plant, which leads to more energy being 
available for the bottoming cycle. A thermal efficiency of 41.5% 
may seem somewhat low, but for an integrated gasification plant 
producing only 120 l<W e , this efficiency is in fact sufficient. Note 
that an integrated gasification combined cycle has an efficiency of 
approximately 40% for a power output of approximately 500 MW. 
Note also that the SOFC fuel utilization is relatively low (0.675), as 
explained earlier. Furthermore, 130.52 l<W t h of domestic hot water 
is produced as a byproduct. 

The cost of the produced electricity is 0.1204 $/kWh, and the 
cost of the produced DHW is 0.0214 $/l<Wh. These costs will change 
if the assumptions given above change. However, the results ob¬ 
tained here provide a fairly good overview of the cost situation. 

The obtained costs are a result of low fuel costs and high 
component and installation costs with respect to the plant size. As 
stated in energy.eu [42], the electricity price for a household with a 
consumption of 7500 kWh/year (±30%) is 0.2562 €/kWh in 
Denmark, and for an industry with a consumption of 2 GWh/year 
(±50%), the price is 0.0982 €/kWh. The results obtained at today’s 
exchange rate from € to $ lead to an electricity cost that is within 
the confines of what is expected when buying it from the grid; 
hence, the system will most likely be competitive for installation in 
places such as hotels and malls. Thus, the system might be cost 
effective when SOFCs and Stirling engines enter the commerciali¬ 
zation phase, and their price will be close to that predicted herein. 
In this study, DHW is considered to be a byproduct, and the cost is 
obtained as 0.0214 $/kWh, which is much lower than the price of 
district heating networks (0.1143 $/kWh as provided by Kobenhavn 
Energi 2012 [43]). Thus, for in-house use, electricity is produced at a 
cost that is nearly competitive when fuel prices are held at the 
assumed level. However, for selling the electricity to the grid, the 
production cost is considered to be high. Compared to other 
renewable energy sources of similar sizes, conventional biomass 
gasifiers at 20—50,000 kW produce energy with costs of approxi¬ 
mately 0.13 $/kWh, while small-scale wind turbines produce en¬ 
ergy at even higher rates [44], 

Here, the total TCI for the installed capacity is 3432.97 $/l<W, 
which is high compared to other renewable energy systems. 


5. Conclusions 

A small-scale integrated gasification SOFC—Stirling CHP plant 
with a net capacity of 120 l<W e is presented and subjected to 
thermodynamic and thermoeconomic investigations. A rather 
modest plant thermal efficiency is found (0.41) after thermoeco¬ 
nomic optimization, partially because the parameter inputs for the 
different components have been chosen from a conservative 
viewpoint and partially because the produced syngas contains a 
large fraction of non-usable compounds, which in turn resulted in a 
lower utilization factor for the SOFC stacks. 

The thermoeconomical analyses showed that by reducing the 
number of stacks from 160 stacks with 74 cells to 120 stacks, the 
cost of the produced electricity decreases at the expense of plant 
thermal efficiency. An electricity production price of 0.1204 $/kWh 
is found with a DHW production cost of 0.0214 $/kWh based on the 
assumption of component cost equations for the future pricing of 
SOFCs and Stirling engines when they emerge into the commer¬ 
cialization phase. These prices are competitive in the Danish mar¬ 
ket for in-house usage but are slightly high if the electricity is sold 
to the grid. 

Neglecting disposal costs, the plant cost is estimated as 3433 $/ 
kW, which is competitive when compared to other environmen¬ 
tally friendly energy systems of similar sizes. 
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Nomenclature 

Ahex heat exchanger area, m 2 

A jj matrix 

C P product cost rate, $/kWh 

C F fuel cost rate, $/l<Wh 

c p specific energy cost, $/l<W 

Dbin binary diffusion coefficient 

D ce n cell diameter, m 

E exergy flow rate, kW 

Epc electricity from fuel cell, V 

£ Nemst Nernst ideal reversible voltage, V 
e specific exergy, kj/kg 

F Faraday’s constant, C/mol 

frf efficiency correction factor 

k thermal conductivity, W/m 2 1< 

g° standard Gibbs free energy, J/mol 

gr Gibbs free energy, J/mol 

h enthalpy, J/kg 

h f enthalpy of formation, J/mol 

f C omp purchase cost of component k 

Lceii cell length, m 

m mass flow rate, kg/s 

hp[ 2 molar reaction rate of H 2 , mol/s 

n e number of electron 

P power, W 

p pressure, bar 

Ph 2 partial pressures for H 2 , bar 

Ph 2 o partial pressures for H 2 O, bar 

Q heat, J/s 

q j interest rate 

T operating temperature, I< 

f thickness, m 

R universal gas constant, J/K mol 

RV reversibility factor 

Up fuel utilization factor, the actual fuel used in the SOFC 
cells compared to full 

V volume, m 3 

i/ an anode porosity 

W work, W 

Y molar fraction 

Z 01 annual contribution to investment, $/kWh 

Z° M annual contribution to operation and maintenance, 

S/kWh 

Greek symbols 
A change/difference 

ATi m logarithmic mean temperature difference, K 

e effectiveness factor 

tj efficiency 

ijpcy polytropic efficiency 

v specific volume, m 3 /kg 

Subscript 

act activation polarization 

an anode 

aux auxiliary 

ca cathode 

cone concentration polarization 


D destroyed 

el electricity 

FC fuel cell 

L loss 

mec mechanical 

Nernst Nernst ideal reversible voltage 

ohm ohmic polarization 

ref reference 

rev reversible 

th thermal 

v voltage 

Abbreviations 

AP anode preheater 

CF1P combined heat and power 

CP cathode air preheater 

DHW domestic hot water 

DNA Dynamic Network Analysis 

EES Engineering Equation Solver 

GAP gasifier air preheater 

CB catalytic burner 

LHV lower heating value 

LMTD logarithmic mean temperature difference 

NTU number of transfer units 

SG steam generator 

PEC component purchase cost 

RH room heater 

SOFC solid oxide fuel cell 

TCI total capital investment 

WH water heater 
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